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Abstract. For understanding carbon erosion and redeposition in nuclear fusion 
devices, it is important to understand the transport and chemical break-up of 
hydrocarbon molecules in edge plasmas, often diagnosed by emission of the CH A 2 A- 
X 2 n Gcro band around 430 nm. The CH A-level can be excited cither by electron- 
impact or by dissociative recombination (D.R.) of hydrocarbon ions. These processes 
were included in the 3D Monte Carlo impurity transport code ERO. A series of 
methane injection experiments was performed in the high-density, low-temperature 
linear plasma generator Pilot-PSI, and simulated emission intensity profiles were 
benchmarked against these experiments. It was confirmed that excitation by D.R. 
dominates at T e < 1.5 eV. The results indicate that the fraction of D.R. events that 
lead to a CH radical in the A-level and consequent photon emission is at least 10%. 
Additionally, quenching of the excited CH radicals by electron impact de-excitation 
was included in the modeling. This quenching is shown to be significant: depending 
on the electron density, it reduces the effective CH emission by a factor of 1.4 at 
n e = 1.3 * 10 20 m~ 3 , to 2.8 at n e = 9.3 * 10 20 m~ 3 . Its inclusion significantly improved 
agreement between experiment and modeling. 
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1. Introduction 

Carbon has a long history as the plasma-facing material (PFM) of choice for the 
walls, limiters, and divertors of tokamaks, including divertor plates of ITER. However, 
upon bombardment with hydrogenic atoms and ions, methane and other hydrocarbon 
molecules are formed, which can easily leave the PFM surface and enter the edge plasma. 
This chemical erosion is very efficient even at low ion impact energies [1] . At the strike- 
points of the ITER divertor, where ion fluxes are maximal, electron temperatures (T e ) 
are expected to be around 1-10 eV, much lower than in current tokamaks. The expected 
electron density (n e ) is approximately 10 20 m -3 , which is much higher than in present 
tokamaks [2], leading to substantial chemical erosion. The hydrocarbons ejected that 
way from the PFM will be transported, ionised and dissociated in the edge plasma. They 
will often be redeposited at other locations, where they may retain unacceptably large 
amounts of radioactive tritium [3] 0]. Both erosion and tritium retention may limit 
the availability of next-generation carbon-clad tokamaks and are therefore of utmost 
importance [5]. 

A common diagnostic for chemical erosion of carbon is the optical emission of the 
CH radical, in particular the emission around 430 nm due to the A 2 A — 7-X 2 n transition 
(the Gero band). Computer simulations are often benchmarked to these measurements, 
using the assumption that the A 2 A level is primarily excited by inelastic electron- 
impact (E.I.) collisions. However, the rate of electron impact excitation decreases 
for lower electron temperatures. Therefore, when electron temperatures are around 
ITER-divertor-relevant (T e < 3 eV) levels, this is no longer necessarily the dominant 
excitation mechanism of the A 2 A level. The exothermic dissociative recombination 
(D.R.) reactions shown in table [l] can be expected to produce, amongst others, excited 
CH radicals in the A-level; this process may dominate at low electron temperatures. As 
previous methane injection experiments [6] in the Pilot-PSI linear plasma generator [7] 
have shown, the absolute Gero band emission efficiency stays almost constant when the 
electron temperature is decreased from 1.2 eV down to 0.1 eV. The cross sections for 
D.R. are only weakly dependent on temperature; therefore those results suggested that 
D.R. is the dominant excitation mechanism. This was also seen in expanding thermal 
plasma [E], in the divertor simulator MAP-II [9], in a magnetically enhanced capacitive 
RF discharge [10] . and in plasma spraying [11]. Understanding of the photon yield from 
D.R. is essential for interpretation of CH plumes (i.e. emission patterns), which would 
give confidence in predictive modeling for ITER plasma-surface interaction and tritium 
retention [TJ1 M, E] • 

The importance of D.R. in those plasma conditions was investigated by comparing 
results of the Monte Carlo impurity transport simulation code ERO [15] with 
experiments in the linear plasma generator Pilot-PSI [33]. The first simulations of 
hydrocarbon transport in Pilot-PSI have been described elsewhere [17] ; the present 
paper is a continuation of those modeling efforts. A number of modifications were made 
to ERO in order to improve its suitability to Pilot-PSI conditions. Several methane 
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Table 1. Dissociative recombination reactions for the methane family, with the labels 
used throughout this paper to refer to them. Also given are their exothermicities, 
taken from [115] 



Label 


Reaction 






Exothermicity (eV) 


CRf D.R. 


e + CH+ - 


-> ch(*; 


) + H + H 2 


3.42 


CH+ D.R. (a) 


e + CH+ - 


* CH(* 


)+H 2 


5.10 


CH+ D.R. (bj 


e + CH+ - 


■» CH + 


H + H 


0.64 


CH+ D.R. 


e + CH+ - 


-> CH(* 


)+H 


6.00 



puffing experiments were analyzed, including two newly presented experiments. It was 
investigated whether the combination of excitation by electron impact and dissociative 
recombination into excited levels can satisfactorily explain the observed photon emission 
intensities in each of these experiments. De-excitation of the CH A-level by electron 
impact was included in the modeling, and shown to be an important process at high 
densities. 

Hydrocarbon dissociation and excitation mechanisms at T e < 2 e V 

At low electron temperatures, rates for ionisation and dissociative excitation by electron- 
impact rapidly go down due to the threshold energies involved. The hydrocarbon 
(methane) molecules entering the plasma are then predominantly dissociated by a chain 
of charge exchange reactions, followed by dissociative recombination [18] . Due to the low 
fraction of hydrocarbon impurities relative to the hydrogen concentration, the reactions 
between two hydrocarbon molecules can be neglected. Charge exchange (CX) is the 
dominant mechanism for ionization of the neutral molecules: 

CH (a0 + H+ -> CH+_ y) + (H, H 2 ) (y < 1). (1) 

There are a number of different dissociative recombination channels that each 
hydrocarbon ion can undergo. For the methane family (CH+), four of these reactions 
can produce CH radicals, relevant for the Gero band emission. They are listed in table [TJ 
Each of these reactions is exothermic, and the CKj~ D.R., CHj" D.R. (a), and CHj D.R. 
reactions are each able to overcome the 2.88 eV electronic excitation threshold of the 
CH radical. Thus, they should be considered as potential populating processes of the 
CH(A) level. The CH^ D.R. (b) reaction can not excite the A-level. 

The probability that each particular reaction actually excites the A 2 A level is not 
known beforehand. There are various other energetically accessible excited states (ie. 
B 2 X~, C 2 £ + ), and the exothermicity can also simply be transferred to kinetic energy of 
the reaction products. A measurement of the population of each of the different excited 
states was not available. The order of magnitude of the excitation probability of A 2 A 
is estimated in section [5] by comparing modeling results to experimental values. 
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2. CH4 puffing experiments 

In various experiments, methane (CH 4 ) was injected into the magnetised hydrogen 
plasma beam of Pilot-PSI. The plasma beams used in this study typically had T e ~ 1 
eV, n e ~ 10 20 m~ 3 at the plasma axis with approximately 10 mm full width half 
maximum (FWHM) for both parameters. The magnetic field was 0.4 T. CH emission 
was measured using an absolutely calibrated CCD camera equipped with a bandpass 
filter (peak transmission at 432.2 nm and a 2.0 nm FWHM Gaussian transmission 
function). Reference light intensities observed without CH puffing were subtracted from 
the observed light. The measured intensity was then multiplied by a factor 2.8 to correct 
for the fact that only part of the CH band passes the filter [20] . This way, line-integrated 
2D profiles of CH emission can be conveniently measured. Deviations from these settings 
are noted with the individual experiments. 

The total emission efficiency is characterised by the photon efficiency 
n = (j)%n x [photons ■ s' 1 ] 
phot y C h 4 [molecules • s _1 ] ' 

where * s the total photon flux per second and Tch 4 is the number of injected 

methane molecules per second. By assuming that the dominant erosion product is CH 4 , 
one can calculate the gross erosion flux in an erosion experiment by multiplying the 
measured CH A — > X photon flux by H p hot- For this reason, the photon efficiency is 
commonly used for calibrating erosion measurements [2D [221 1231 [21] • Under the "corona 
assumption" (a homogeneous low-density plasma), H p hot is determined only by the rate 
coefficients of hydrocarbon dissociation and of CH excitation. In such cases it can be 
obtained by 0-dimensional calculation, such as the HYDKIN toolbox [25J. But in reality, 
carbon redeposition and reflection as well as temperature gradients and transport losses 
in the plasma do influence the emission significantly. This makes detailed transport 
modeling necessary for CH emission and erosion quantification in tokamaks. 

3. Simulation details 

The present study is a continuation of previous modeling efforts in Pilot-PSI [T7] . 
Methane test particles leave a puffing hole (0 0.6 mm) in the target at a thermal 
energy distribution corresponding to a source temperature of 700 °C (approximately the 
target temperature), with a cosine angular distribution. Quickly after injection, they 
get ionised and dissociated predominantly through charge exchange and dissociative 
recombination. Reaction rates for methane breakdown come from [19J. The modeling 
used a static background plasma, with plasma parameters from experimental Thomson 
scattering measurements. 

Two photon emission channels are taken into account. Cross sections for CH 
excitation by electron impact were taken from [26j. Furthermore, every dissociative 
recombination event capable of producing excited CH was counted (keeping in mind 
that only a fraction of these will actually produce excited CH). In both cases, light 
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emission is assumed to occur at the same point in space as the excitation, due to the 
short lifetime of the CH A-level. Due to the high electron densities, it appeared necessary 
to also take de- excitation by inelastic electron impact collisions and by dissociation into 
account. 

Though this study focuses on the quantitative description of photons from 
dissociative recombination, several improvements have been made in comparison with 
previous modeling. They will be outlined below. The effects of these modeling 
assumptions on simulation results are discussed in section [6] along with an error analysis. 

3.1. De-excitation of CH (A 2 A) 

The radiative lifetime of the CH A 2 A level is relatively short (0.53 /is [27]). Therefore, 
in low-density plasmas it is often assumed that excitation of a CH molecule is always 
followed by radiative decay. However, in high-density plasmas, one should not neglect 
the importance of de-excitation of CH by electron impact and by dissociation. The rate 
of electron de-excitation can be obtained from the excitation rate by using the principle 
of detailed balance [28J: 



with kx~>A an d k^x the excitation and de-excitation rates, g the statistical weight of 
the given levels (where for CH, g& = gx), and AE the excitation energy. At higher 
electron densities, this de-excitation causes quenching of the photon emission. 

There is more than one mechanism that can quench the CH A— > X photon emission. 
Quenching also occurs if the rate of ionisation or dissociation is large enough to be 
comparable with the inverse radiative lifetime of CH. This possibility was taken into 
account, but found to be of relatively minor importance, as the rate of dissociation 
is much lower than that of electron impact de-excitation. The rate coefficient of the 
primary dissociation channel of CH (charge exchange with H + ) is 1.29 ■ 10~ 15 m~ 3 at 
Tj = 1 eV, resulting in an average lifetime of 8 (is, which is one order of magnitude 
above the radiative lifetime. It is also possible that electron impact causes transfer 
from the A 2 A level to other excited states, such as the B 2 £~ or the C 2 S + levels. In 
particular, it is known that the A 2 A, v — and the B 2 £~, v — 1 are near-degenerate 
and collisional interconversion between these states can proceed efficiently [29]. Since 
the population of the B 2 S~, v = 1 level is typically much lower than that of the A 2 A, 
v = level, such collisional coupling should cause net reduction of the A 2 A, v = 
population. A complete collisional-radiative model for the CH radical would take this 
into account, but unfortunately no reaction rates are available for the electron-impact 
collisional interconversion between these levels. Therefore the quenching rate used herein 
should be seen as a lower limit of the actual quenching rate. 

The A-level is depleted by both spontaneous emission (occurring at a rate of the 
Einstein A-coefficient Aa^x) and by de-excitation (occuring at a rate n e ■ kA^x)- 




(3) 
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Therefore the effective photon emission rate is reduced by the same factor: 

A A - 



f. 



quench 



(T e ,n e 



rX 



n, 



(4) 



kA^x + A A ^x 

such that fquench = 1 indicates no quenching, and f qU ench < 1 indicates that quenching 
reduces the effective photon emission. The quenching factor is plotted in figure [T] 
Electron densities at the strike point of the ITER divertor are expected to be above 
10 20 m~ 3 . In such conditions f qU ench is significantly below 1, such that CH quenching by 
de-excitation should not be neglected. 
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Figure 1. Quenching factors for different electron temperatures. 



3.2. Pre-sheath treatment 

Plasma ions are accelerated towards the sheath edge in the so-called pre-sheath region. 
The sheath itself is very thin (~ 10 /im), but the pre-sheath length is on the order of 
a cm. Typically, the decay length of the CH emission intensity is in the same order of 
magnitude. Therefore it is necessary to consider the acceleration and associated density 
drop in the pre-sheath. Due to the low electron temperatures in Pilot-PSI, ionisation 
inside the pre-sheath is negligible. Pilot-PSI has a flowing plasma, and collisions with 
neutral hydrogen are significant; this strongly influences the pre-sheath. In Pilot-PSI, 
the pre-sheath is formed by momentum loss due to collisions between plasma ions and 
the neutral hydrogen before the target. Such a "collisional pre-sheath" has different 
properties than the ionizing pre-sheath which is found in plasmas with higher electron 
temperatures [50] . 

The thickness of the pre-sheath is determined by the mean free path X n i with which 
incoming plasma ions lose their momentum to the neutrals before the target. Typically, 
the main momentum loss mechanism is charge exchange with H atoms reflected from 
the surface. Charge exchange with molecular hydrogen and elastic collisions are only of 
secondary importance. A typical value for X ni of 7 mm has been found by fitting model 
and experiments [31]. The corresponding H density nn can be found using the cross 
section o~cx for the dominating charge exchange reaction (H + H + — > H + + H), which 
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is 6 ■ 1CT 19 m 2 [32J. From that, one finds n H = l/(A ni * o~cx) ~ 2 ■ 10 20 m~ 3 . The ion 
flow velocity far away from the target was taken from measurements [33J to be 4000 
m s . The density drop and velocity profile corresponding to those parameters were 
implemented in these simulations. The resulting density drop between the Thomson 
scattering laser and the target is a factor of 3.4. This number is somewhat greater than 
the typically assumed pre-sheath density drop of a factor 2, because the friction between 
ions and neutrals is important in Pilot-PSI. 

3.3. Other modeling improvements 

The Pilot-PSI plasma flow normally rotates around the plasma axis. Though maximum 
rotation velocities of over 10 km s _1 were measured close to the plasma source |34J, 
rotation velocities measured close to the target were significantly lower [31]. In this 
work, the maximum rotation velocity was set at 3000 m s -1 at r = 7.0 mm away from 
the plasma axis; this value was taken from the experiment best matching conditions 
in this paper. The electric field that causes this plasma rotation was also taken into 
account. No big qualitative changes were found due to this rotation. 

In Pilot-PSI the ion density is of the same order as the neutral density [35] . 
Therefore Coulomb friction is the dominant force affecting the test particles. For 
that reason collisions between the test particles and neutral hydrogen were neglected, 
speeding up the simulations. Because of the strong re-erosion caused by the high 
hydrogen flux to the target, carbon redeposition was neglected. All carbon molecules 
arriving at the target are assumed to reflect as CH4. 

4. Results 

4-1. Puffing into the side of the plasma beam 

Methane was injected perpendicularly through a nozzle into the side of the plasma 
beam at a rate of 3 seem (1.34 • 10 18 molecules s -1 ). The injection nozzle (0.6 mm 
diameter) was located 25 mm away from the target and 17 mm away from the center 
of the beam. There, the axial gradients of the plasma parameters are small. This setup 
has two advantages. First, carbon recycling at the target is not so important here, due 
to the distance between the nozzle and the target. Second, the Thomson scattering 
laser that measures electron temperature and density is located at the same distance 
from the target as the CH4 injection nozzle, thus enabling measurements of the plasma 
parameters at the location of injection. 

Both the measurements and the simulations were performed with two different 
electron temperatures, which we will call the "high T e " and "low T e " cases. Electron 
temperature and density, as measured by Thomson scattering, are shown in figure [2j 

Measured and simulated plumes are shown in figure [3j We do not yet know which 
fraction of the D.R. events will produce a CH radical in the A 2 A level (see section [5]). 
Therefore, the D.R. "photon count" is an upper limit for the actual number of photons. 
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R (mm) R (mm) 

Figure 2. Thomson scattering measurements for the side puffing experiment, and the 
Gaussian fits that were used for simulation input. 

It is clear from these results that the dissociative recombination produces a qualitatively 
different emission plume than the electron-impact (E.I.) excitation. The total photon 
yield from E.I. excitation is also far less than that from D.R., even in the experiment 
with higher T e . The location of the simulated D.R. plume is in good agreement with 
the experimentally observed plume size, indicating that in both conditions electron- 
impact excitation is not the dominant process. Note that this does not imply that the 
predicted electron-impact excitation is in mismatch with the experiment; it is simply 
not the dominant process under these conditions. 

Another interesting observation is that, similar to earlier methane injection 
experiments from a different geometry [6], the observed photon emission coefficients 
were almost independent of electron temperature. 

Table 2. Effective photon emission coefficients Uphot as measured in the side-puffing 
experiment and modeled values from different processes. 



Kphot (low-T e case) U phot (high-T e case) 



Experiment 


0.038 


0.033 


E.I. excitation 


<0.00I 


0.009 


CH+ D.R. 


0.003 


0.019 


CH+ D.R. 


0.012 


0.032 


CH+ D.R. 


0.025 


0.101 



The total effective photon emission coefficients Uphot (@) are given in table [2] Of 
the D.R. channels, the CfLj D.R. (e + CH| -»■ CH(*) + H + H 2 ) occurs by far most 
frequently. This does not necessarily imply that the reaction will be dominant in actually 
exciting the A-level, since the probability of yielding an excited radical is not given, as 
will be discussed in section |5j However, the sum of all other emission channels is not 
enough to explain the observed photon emission yields. That means that the fraction 
of CH4 D.R. that yields an excited radical should be significant in order to explain the 
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Figure 3. Experimental (left) and simulated (right) emission profiles during the side- 
puffing experiment. In the high-T e case, the modeled number of photons from both 
electron-impact (E.I.) excitation and dissociative recombination (D.R.) are shown. E.I. 
excitation is negligible in the low-T e case. The plots of D.R. events show the sum of 
all relevant D.R. channels. The small light plume close to the injection nozzle in the 
low T e experiment is due to reflection on the injection nozzle. 
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Illustration of the experiments of section 4^2 and section 4.3 



experimental H p hot at low T e . 

4-2. Puffing into the center of the plasma beam 

In the following two series of experiments, methane was injected directly into the center 
of the plasma beam, through a hole in the target. This geometry is illustrated in figure [4j 
Profiles of the (line-integrated) CH emission in the center of the plasma beam are 
shown in figure [5] together with simulation results. The profile is taken along the dashed 
line in figure |4| On a logarithmic scale, the light emission quickly drops after several 
millimeters. From the decay length of the various excitation channels, one can see that 
the dissociative recombination of CH3" gives a good match with experiment. Also the 
DR of CHj will match experiment, if weighted by an appropriate excitation probability. 
CHg" DR, by itself, is not enough to explain experimental photon emission. 




o 

eu 1E14 



2 4 6 

Distance from target (mm) 



Figure 5. Photon emission profiles in the center of the plasma at T e =L25 eV. 
Experimental values are from [6]. 



We would now like to establish the T e -dependence of the various photon emission 
processes. The absolute CH A-X photon emission was measured during scans of T e . 
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These experiments are described in detail elsewhere: "experiment A" was performed 
using a spectrometer [6] , and "experiment B" was performed using a CCD camera [36J . 
Both the spectrometer and the CCD camera were absolutely calibrated. Therefore, 
both experimental results can be compared with each other and with simulation results. 
Results are shown in figure [6j The photon yield of each of the D.R. processes (table 
[TJ, as well as from electron impact excitation is shown separately. It should be noted 
that the experimental point at 2.4 eV had a higher magnetic field (0.8 T, rather than 
0.4 T in the other points); this led to a relatively high electron density. That might 
explain the low photon yield in that experiment. Nevertheless, the magnitude of this 
discrepancy is unexpected, and, for the time being, unexplained. 

Below 1.5 eV, E.I. excitation drops steeply due to the threshold character of that 
process. In that regime, the dissociative recombination channels are clearly necessary 
to explain the photon yield. Like in section 4.1 CHj dissociative recombination is 
the most frequently occurring D.R. channel. By just counting CHj D.R. events, we 
overestimate the experimental H p hot- This suggests that only a fraction of these D.R. 
events actually produces a photon. Due to the discrepancy between experimental results 
at T e > 2 eV, it is not yet possible to give a precise temperature threshold at which 
D.R. becomes dominant. We can, however, conclude that certainly below T e = 1.5 
eV, D.R. is necessary (and sufficient) to explain the intensity of CH emission in these 
experiments. 
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Figure 6. Tl p hot from experiment and from simulation of the various CH excitation 
channels, as a function of T e . Details can be found in the main text. 



4-3. Electron density scan, importance of CH de-excitation 

A scan of the electron density was performed both experimentally and in modeling. 
Such a scan can shed light on the role of de-excitation of the CH emission (see section 
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3.1). This quenching increases as the electron density goes up, as shown in figure [T] 
In these experiments, methane was puffed through the target into the center of the 
plasma. As in [6], the photon yield was measured with a spectrometer; however, it 
has not been absolutely calibrated. The magnetic field was set to 0.8 T to reach high 
electron densities. T e varied between 0.5 and 0.9 eV; each individual simulation used 
the T e and n e that were measured in the corresponding experimental point. 

In one simulation series the effect of de-excitation was taken into account, in another 
series it was neglected. Both were then compared with experiment. Because T e < 1 
eV, the contribution from E.I. excitation is relatively small. Hence, only the sum of all 
dissociative recombination events was compared with experimental values. Results are 
shown in fig. [7| Because the spectrometer was not absolutely calibrated, all results are 
normalised to 1 at the lowest density point. 

It should be noted that the density was measured by a Thomson scattering 



diagnostic at 25 mm before the target. As described in section |3T2] the density drops in 
the final cm before the target, due to friction with neutrals in the pre-sheath. 
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Figure 7. Normalised photon yield from experiment, and D.R. photons from ERO 
simulation at various densities. D.R. photons are summed from all relevant channels 
in table [TJ Densities on the x-axis are those measured at the location of the Thomson 
scattering diagnostic. 



At the lowest electron density (n e = 1.3 * 10 20 m~ 3 ), the calculated emission 
rate is already reduced by a factor 1.4 due to the quenching (though this is hidden 
by the normalisation). At the highest electron density (n e = 9.3 * 10 20 m -3 ), the 
reduction by quenching is a factor 2.8. Experimentally, it was observed that the 
photon yield goes down as the electron density increases, and this was attributed to de- 
excitation. However, the simulations without the inclusion of de-excitation give exactly 
the opposite result: the photon yield goes up with increasing density. The explanation 
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for this modeling result is that with a higher electron density, hydrocarbon molecules 
are dissociated more rapidly, significantly reducing the loss fraction of hydrocarbon 
molecules that leave the plasma before they are dissociated to CH. The simulation 
series that includes CH de-excitation is considerably closer to experiment than the series 
without de-excitation. 

Even though the trend of a decreasing photon yield which increasing density is not 
yet fully reproduced, these results do underline the significance of de-excitation of CH 
in high-density plasmas. Since collisional interconversion between the A 2 A level and 
other excited levels could not be taken into account, the quenching rate used in these 
calculations is a lower limit of the actual quenching rate (see section 3.1). That is a 
possible explanation for the discrepancy between the experiment and simulation. 



5. CH A-level excitation probabilities of the different dissociative 
recombination channels 

In all of the simulations, CH4" D.R. is the most frequently occurring D.R. channel. This 
can be explained by examining the reaction database that is used by ERO |19j . From 
these reaction rates it follows that many of the methane molecules entering the plasma 
are quickly ionised to CHj by charge exchange with H + . CH4" then quickly recombines 
to a variety of products, including (25 %) CH. The other D.R. channels that can excite 
the CH radical are CH^ D.R. and CH^ D.R. 

From these results, it is not trivial to estimate the fraction of each D.R. channel that 
excites the CH level. The electron temperature scan in figure [6] suggests that in the case 
of CH^J~ D.R. this fraction should be on the order of 10% to match the experimentally 
observed photon flux. However, looking at table [2j one estimates a somewhat larger 
value for this fraction. Either way, it does seem that excitation of the CH A 2 A level is 
relatively common during the dissociative recombination of CH+ ions. Both experiments 
suggest that these probabilities are at least of the order of 10%. Judging from figure [6j 
it seems that the effect of E.I. excitation does become significant starting from electron 
temperatures around 1.5-2.5 eV. Thus, even though D.R. excitation of CH probably 
does generate some photons in the CH4 puffs in detached plasma at JET (T e ~ 2 eV) 
[37J and DIII-D (T e = 2.0 eV) [38J, it might be a marginal process. Further research is 
needed to conclusively answer this question. 

Unfortunately, the branching ratio of the dissociative recombination of CH^J~ into 
various hydrocarbons CHj, has never been directly measured. Instead, these branching 
ratios are extrapolated from measurements on D.R. of CH^~ [39] and CHj^~ jlQ]. Because 
of the apparently important contribution of this channel at low T e , experimental 
confirmation of this branching ratio would be beneficial to modeling efforts. 

The fact that our simulations consistently give the D.R. of CH^ to be the most 
frequently occuring process that can excite CH complicates the benchmarking of 
spectroscopy during erosion experiments with methane injection experiments at low 
T e . In models of chemical erosion, CH 3 is often assumed to be an important eroded 
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hydrocarbon species jH]. That species does not readily form CH4" in the plasma. This 
means that the dominant excitation mechanism may be different in erosion experiments 
and in injection experiments, complicating analysis of erosion experiments. 

6. Discussion 

Aside from the general measurement error in the underlying reaction rate databases, 
there are a number of uncertainties specific to these simulations; the major ones will be 
summarised briefly here. 



The properties of the pre-sheath (discussed briefly in section 3.2) have a significant 



influence on the results of ERO calculation, which is greater than the effect of other 



assumptions on plasma parameters (section 3.3). The drop in electron density close 
to the target causes the impurity penetration depth to increase to a greater value 
than the experimental observation. However, varying the pre-sheath parameters 
within a realistic range does not have such a great influence on the integrated 
photon emission yield; this remains within a factor two throughout a parameter 
scan. For that reason, only integrated photon yields are compared here. The only 
exception is the puff into the side of the beam, which is not so strongly affected 
by the pre-sheath. The exact pre-sheath density profile and plasma acceleration 
remains an interesting subject for future study. 

With puff-rates of up to 3 seem, the injection flux density right behind the target 
can be comparable to the hydrogen flux. For instance, when puffing into the center 
of the plasma through a 0.6 mm nozzle, the hydrocarbon flux density right behind 
the nozzle is 1 • 10 24 m~ 2 s -1 , which is about the same as the incoming hydrogen ion 
flux. Furthermore, cauliflower-shaped carbon microparticles have been found on 
the target both in carbon erosion and puffing experiments in Pilot-PSI (12} l4~3l l4~4"] . 
If these are formed in the plasma, it means there has to be formation of higher 
hydrocarbons (C^H^, x > 1). Therefore one should investigate whether the test 
particle assumption (i.e. no collisions between impurity molecules, and no influence 
of the impurity on the plasma) is still valid. If not, the influence of the hydrocarbon 
injection on local plasma parameters would have to be taken into account, as was 
done earlier for injection experiments in TEXTOR [33]. 

A scan of the injection rate from 0.6 to 3 seem has been performed [3S]- There, no 
obvious disturbance of the plasma, change in the CH penetration depth, or non- 
linearity in the total photon emission yield was observed. Thus, it seems that the 
central assumption of ERO is reasonable. The explanation for that may be that 
hydrocarbon molecules quickly spread radially through the plasma beam, so that 
their density is only a significant fraction of the plasma density in a very small area 
behind the injection nozzle. 

Unlike previous calculations [17], all hydrocarbon radicals and molecular ions 
arriving at the target were assumed to be reflected. The reason for this is the 
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high hydrogen flux in Pilot-PSI, which causes strong re-erosion of the carbon 
deposition; normally, only a small fraction of the injected carbon is retained on 
the metallic target. Reference simulations were also performed with the opposite 
extreme assumption: all carbon molecules arriving at the target are permanently 
deposited. That extreme assumption typically reduced the calculated total photon 
yield by approximately a factor 2 compared to the simulations in this paper. The 
penetration depth of carbon in the plasma is not significantly affected. We conclude 
that the uncertainties about sticking do not affect qualitative results. 

7. Conclusions 

Many updates have been made in simulations of methane puffing experiments in Pilot- 
PSI. Most importantly, the contribution of dissociative recombination to the production 
of excited CH has been included in the modeling. By comparing simulations to a series 
of both new and existing CH 4 injection experiments, it was then confirmed that this 
process dominates at electron temperatures below 1.5 eV. These results suggest that 
the fraction of D.R. events that actually produce an excited CH radical is at least of 
the order of 10%. 

De-excitation of CH* by electrons was included in the modeling, and an 
experimental scan of electron density was performed. It was demonstrated that at higher 
electron densities, such as expected in the ITER divertor, de-excitation certainly cannot 
be neglected. Inclusion of photon emission quenching by electron-impact de-excitation 
in the modeling greatly improved the agreement between modeling and experiment in 
a scan of n e . 

In present tokamaks, electron temperatures near the plasma-facing surface are 
typically much greater than 3 eV. Under such conditions, electron-impact processes 
are very efficient in the dissociation of hydrocarbon ions. In simulations of such 
plasmas, good agreement is typically obtained between measured and predicted photon 
emission [16], and the photon production by dissociative recombination can be neglected. 
However, this is clearly not the case for low-temperature, ITER divertor-relevant 
plasmas. 

These observations shed new light on the interpretation of CH photon emission 
measurements. Finally, they could also be of relevance for low-temperature hydrocarbon 
laboratory plasmas experiments, which also often employ CH spectroscopy. 
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